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Using mode-coupling theory, we derive a constitutive equation for the nonlinear rheology of dense 
colloidal suspensions under arbitrary time-dependent homogeneous flow. Generalizing previous re- 
sults for simple shear, this allows the full tensorial structure of the theory to be identified. Macro- 
scopic deformation measures, such as the Cauchy-Green tensors, thereby emerge. So does a direct 
relation between the stress and the distorted microstructure, illuminating the interplay of slow struc- 
tural relaxation and arbitrary imposed flow. We present flow curves for steady planar and uniaxial 
elongation and compare these to simple shear. The resulting non-linear Trouton ratios point to a 
tensorially nontrivial dynamic yield condition for colloidal glasses. 



PACS numbers: 82.70.Dd, 64.70.Pf, 83.60.Df, 83.10.Gr 



The imposition of flow can drive a fluid far from equi- 
librium. Due to the occurancc of long relaxation times, 
this effect is ubiquitous in complex fluids (colloids, poly- 
mers, etc.), whose rheology is of significant technologi- 
cal interest, and also represents an important challenge 
in nonequilibrium statistical physics. Continuum ap- 
proaches have provided important insights, using sym- 
metry and other principles to construct or constrain phc- 
nomenological constitutive relations. While the constitu- 
tive equations of Newtonian fluids and Hookian solids are 
derivable from fundamental starting points (the theory 
of linear response based on Onsager's regression hypoth- 
esis), there has been less progress with their nonlinear 
generalizations for viscoelastic fluids, plastic solids and 
other strongly deforming soft materials. 

A central aim of theoretical rheology is thus to derive 
from the underlying microscopic interactions the con- 
stitutive equations that relate the stress tensor to the 
macroscopic deformation history of a material. For en- 
tangled polymer melts, the constitutive equation of Doi 
and Edwards 0| has enjoyed considerable success. An 
analogously general microscopic constitutive equation for 
colloidal dispersions remains conspicuously lacking 
Even the simplest hard-sphere colloids in concentrated 
suspension exhibit a broad range of viscoelastic behav- 
ior; alongside to flow-thinning [3| and thickening [i[ , slow 
structural relaxation leads to glasses showing solid-like 
response, strain- hardening/softening, and plastic flow Q. 

But, while the linear viscoelastic spectra of colloidal 
suspensions are fairly well understood [f|, only recently 
has progress been made in nonlinear flow predictions for 
simple shear 0, H|. Shear represents a relatively weak 
flow in which material lines grow linearly with time, 
while in elongational flows such growth is exponential, 
creating much more severe deformations of material ele- 
ments. Thus a description capable of handling arbitrary 
deformation histories is highly desirable. In the contin- 
uum approaches, invariance arguments strongly restrict 
the deformation measures that can appear. A good mi- 
croscopic theory (e.g., 0|) should implicitly respect such 
invariances so that a tensorially admissible constitutive 
structure emerges from a first-principles starting point. 



In this Letter we develop such a theory for dense sus- 
pensions of spherical colloidal particles under imposed 
time-dependent flow. Our approach, which considerably 
generalizes Q, is nonlinear in the velocity gradient 
tensor and, in the absence of flow, features a transition 
to a glassy solid as a function of the thermodynamic 
control parameters, allowing the delicate interaction be- 
tween glass formation and (time-dependent) external de- 
formation to be investigated. Our treatment assumes in- 
compressible homogeneous flow, neglecting fluctuations 
of the solvent velocity field and thus hydrodynamic inter- 
actions 0. These assumptions may break down at high 
flow rates and/or densities where lubrication forces drive 
cluster formation and shear thickening Q. The nonlin- 
ear theory is formulated by application of mode-coupling 
type approximations to exact generalized Green-Kubo re- 
lations. We thus consider it to be of 'dynamical mean 
field type'. 

We consider a system of N spherical Brownian parti- 
cles dispersed in a solvent with a specified velocity profile 
v(r,i) = n(t) ■ r. The time-dependent velocity gradient 
tensor n(t) is assumed spatially constant, thus excluding 
the inhomogeneous flows which occur in shear-banded 
and shear-localized states. Incompressibility implies that 
Tr n(t) = 0. The distribution function evolves according 
to a Smoluchowski equation dt^(t) = — Y1l% ii (where 
ji is the probability current of particle i) [i|, 0]: 



n(t) 



Wt) = n(t)*(t) 



(1) 



F, is the force acting on particle i due to the other par- 
ticles. The thermal energy fc^T and infinite-dilution dif- 
fusion coefficient Dq are set equal to unity. Equilibrium 
and non-equilibrium solutions to ([1]) are distinguished 
by the existence of a finite probability current. Its exis- 
tence rules out the possibility of a distribution function 
of Gibbs-Boltzmann form. 

When constructing constitutive equations based on 
a given microscopic dynamics it is important to first 
identify the relevant affine deformation measures. We 
achieve this by considering the translational invariance 
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of the two-time correlation functions. The first step is 
to show that a translationally invariant initial distribu- 
tion function can lead to a translationally invariant, but 
anisotropic, distribution function ^(i), despite the fact 
that the Smoluchowski operator is itself not translation- 
ally invariant. This follows from the formal solution 



*(r,i) = e^ s0(r ' s) vMr), 



(2) 



where we assume an equilibrium Gibbs-Boltzmann distri- 
bution ^(T) at t = OjT = {i"i, ■ ■ • , r^} and e + is a time- 
ordered exponential [T(|. We now consider shifting the 
particle coordinates by a constant vector r' { = rj+a. This 
leads to Q(T',t) = 0(r, 0i-/c(*)-a = fl(T,t)+A(t) 

and hence the shifted distribution function 



nr / ,t) = e^ ds(n(r ' s)+A{s H e (T), 



(3) 



where ^ e (T') — ^ e (r). The time-ordered exponential 
can be rewritten using an operator identity 



*(r',i) =ex P+ 



f*dsA(s) 



*e(r). 



(4) 



Using A(i)ty e (r) = and applying the commutator 
[fi(T, s),A(s')} = -a • k t (s) ■ k t ( s >) ■ £. di yields 



*(r / ,t) = *(r,t) > 



(5) 



as desired. We now use this result to study the invari- 
ance properties of the correlation functions. The corre- 
lation of two wavevector-dependent fluctuations is given 

b y C U9 k (MO = ( a /q(*)^ffk(*')>K(t'). wncrc (->K(t') indi- 
cates an average over the distribution (|2|). Translating 

the particles by a constant vector a and using ([5]) yields 



-i(q(t,«') 



(6) 



Translational invariance requires that the correlation 
function is unaffected by the shift with a. This leads 
to the requirement that a fluctuation at wavevector k = 
q(t, t') at time t' is correlated with a fluctuation with 
wavevector q at time i as a result of the affine solvent 
flow, where 

q(M') = q-: />K(s) - (7) 
For shear ||7J) recovers the familiar shear-advection Q, 
while for an extensional flow ([7]) describes the exponen- 
tial deformation of material lines. The tensorial expo- 
nential function in ||7J) may appear unfamiliar, but is 
simply the inverse of the deformation gradient tensor 
F(t,t') = dr(t)/dr(t'), a standard quantity in elasticity 
theory used to connect initial (f') and final (t) coordinates 
following homogeneous deformation, r(t) = F(t, t')-r(t'). 
The deformation gradient is related to the velocity gradi- 
ent tensor via d t F(t,t') = K,(t)F(t,t'). We can thus de- 
fine the forward and reverse-advected wavevectors q(£, t') 
and q(i,t') using the deformation gradient 



The magnitudes of the advected wavevectors are thus 
related to the left and right Cauchy-Green tensors, 
given by B(t,t') = F(t, t')F T (t, i'land C" 1 ^,*') = 
F- 1 (t,t')(F- 1 (t,f)) T , respectively [ll|: 

g 2 (i,t')=q-B(t,t')-q , q 2 (t 7 t') = q-C- 1 (t,t')-< l . (9) 

By considering the translational invariance of solutions 
to (TT]) we have thus identified the appropriate deforma- 
tion measures for describing the affine deformation in our 
Brownian system. This enables us to formulate a theory 
obe ying the principle of material objectivity; see below 
and [lOj . Non-affine particle rearrangements arising from 
excluded volume constraints at finite colloid density are 
approximately treated by the mode-coupling approxima- 
tions to be introduced below. 

Integration over the flow history yields an alterna- 
tive solution of (fT]) which is formally equivalent to ([2|) 
but more suitable for approximation. This integration- 
through-transients approach yields @, [1] 



) *(t) = * e +f dt' * e Tr{Av(t')<7} 

J — OC 



f t *,<fant( a ) 



(10) 



where a a p = — F" r i i s t ne potential part of the mi- 
croscopic stress tensor, and the adjoint Smoluchowski op- 



erator is Cv(t) — 



Yi ■ K T {t)) d,. We take 



the system to be in quiescent equilibrium in the infinite 
past and thus neglect possible non-crgodicity in the ini- 
tial state and related ageing phenomena. The validity of 
this assumption will be dependent upon both the specific 
system and flow history under consideration 0] . By using 
(fTU|) to calculate the average of p* k p k /N and, in the spirit 
of mode-coupling theory, approximating by projecting 
onto density fluctuations Il2| we obtain an expression 
for the distorted structure factor Sk(t; k) = {p%p k ) K {t) 



S k (f,K) = S k - f dt' 

J — OC 



dS, 



k(t,t' 



dt' 



dt' 



, dSk f dq dS, 



q(t.t') OQ 



dn I 16tt 3 dt' S 2 n V q dn q(M ' )V 



k(t,t') 

dS, 



o(M') 



(11) 



where n is the number density and S k the equilibrium 
structure factor used to proxy the colloidal interactions 
fToj | . Eq. lfTT]) describes the flow induced microstructural 
distortion which becomes appreciable when the flow in- 
terferes with slow cooperative structural relaxation [l3l ]. 
Flow enters the description via the advected wavevector, 
both explicitly through Syt,t') an d implicitly through its 
effect on the transient correlator. The second term in (fTTj) 
is anisotropic, whereas the third term is purely isotropic. 

Eq. pop also can be employed to directly approximate 
the stress tensor. Projection operator approximation of 
the resulting generalized Green-Kubo relation provides 
an explicit approximation for <r(t). This we find to be 
related directly to Sk(t; k) by [14[ 



fdkp kk t so / 

q(f , t') = q • F- 1 (t, 0, q(«, f ) = q ■ F(t, t 1 ). (8) °W = J 16^ ~k ° k SS ^ K) ' 



(12) 
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where is the direct correlation function, c' k = dck/dk 7 
5Sk(t] k) = Sk(t; K,)—Sk and II is the equilibrium osmotic 
pressure. (Note that (fT2j) also provides the pressure un- 
der flow). For shear flow cr xy {t) from (fT2"|) coincides with a 
result of Fredrickson and Larson [ll| for sheared copoly- 
mers, reflecting the Gaussian statistics underlying both 
their approach and our own. Stresses are thus connected 
to microstructural distortions, which build up over time 
via the affine stretching of density fluctuations competing 
with structural rearrangements encoded in ^(i, t'). 

In order to close our theory we require the transient 
correlator <J>k(M') = (p£exp_ 

describing the decay under flow of thermal density fluc- 
tuations, where (•) denotes an equilibrium average. The 
appearance of the advected wavevector ensures that triv- 
ial decorrelation due to affine motion alone is removed. 
Time-dependent projection operator manipulations yield 
exact results for the equation of motion of the transient 
density correlator containing a generalized friction ker- 
nel which is amenable to mode coupling approximations 
0, HI . Applying mode coupling type approximations to 
these formal results yields the equation of motion 



d_ 

dt 



* q (Mo) + r q (t,i ) *q(*,fo) 



(13) 



jf di / m q (t,i , ,to)^* q (t' J to)) =° 



where r q (£, to) = Q {t,to)/Sq(t t t )- The friction kernel 
m q (t, s,to) is the autocorrelation function of fluctuating 
stresses, which in mode coupling are connected to struc- 
tural relaxation as described by the density correlator. 
In the present approximation it is given by 



mJt,t',t o y 



dk 



5' ( 3(t,to) 5 'fe(t',to) 5, !5(t',«o) 



(14) 



16tt37 ^(t',£o)g 2 (t,£ ) 
x Wt', *o) Wi, *b)*£(f O*p(t',to)(*. *0. 

Vqk P (t, t ) = q(t, t ) • (k(t, *o)cs( t)to ) + p(t, t )c p ( t)to) ), (15) 

where p = q k. Equations (fTTj) - (fl"5l) form a closed 
constitutive theory for the microstructure and stress re- 
sponse of Brownian particles under external flow, requir- 
ing only Sk and n(t) as input 

We next consider our theory from the st andp oint of 
continuum 'rational mechanics' approaches [ll| . show- 
ing that (in common with [U) it complies with their in- 
variancc requirements but avoids their oversimplifying 
assumptions. Such approaches often express <x(t) as a 
functional of a suitable deformation measure, weighted 
by a fading memory; an example is the integral form 
of the upper-convected Maxwell equation (Lodge equa- 
tion) with the memory taken as a known sum of decaying 
exponentials fl6j . In the present treatment, the mem- 
ory is instead given by $^(t,t') which for general flows 
is both anisotropic and a function of two times rather 
than a simple time difference. Both aspects are physi- 
cally reasonable extensions of the assumption of a fixed 



isotropic decay. In the absence of flow, 4>k(t) exhibits 
non-exponential decay and, in the glass, arrests to a fi- 
nite plateau value leading to solid-like response [l2j . 

An important symmetry consideration is the principle 
of material objectivity (or frame indifference) [H, QJ, E3 ■ 
This asserts that the relationship between er and k should 
be invariant with respect to time-dependent rotation of 
cither the material sample or the observer. Although 
this is only an approximate symmetry, based on the ne- 
glect of inertial effects at the microscopic level [l?} , many 
soft materials display this invariance to a good level of 
approximation. The overdamped Smoluchowski dynam- 
ics ([T]) underlying our treatment excludes inertial effects 
from the outset, so that our set of equations (fITj) -(fl"5 |) are 
material objective, so long as this is preserved by our ap- 
proximations. This can be explicitly confirmed by using 
(HJ and ([9]) to eliminate the advected wavevectors from 
([TlTl- fTS} in favour of the deformation tensors [Toj . 

Next we address the comparison between elongational 
(planar or uniaxial) flow, and shear. Starting with the 
limit of small strain, we obtain the linear response result 



\t) = dt' 



dk 

mT 3 



{(k-/c(t')-k)kk} 



kS k 



(16) 



where linearity enables us to introduce the symmetrized 
tensor R,(t) = («(t) + K T (t))/2, and $fc(t) is the qui- 
escent correlator. As all anisotropy in (|16p is con- 
tained within the factor {•}, the angular integrals may 
be easily evaluated for a given n(t). It follows that the 
planar-cxtensional and shear viscosities are related by 
Vc/Vs = 4, in compliance with one of Trouton's rules. 
In the glass phase the transient correlator does not re- 
lax to zero. Partial integration of (|16|) . followed by tak- 
ing the small-strain limit, leads to <x(i) = 2Ge(t), where 
e(t) is the infinitesimal strain tensor (Tre(t) = 0), G = 
(? f^dk k 4 (c' k cj)k (oo)) 2 / (607T 2 ) is a known mode-coupling 
approximation for the linear shear modulus 0, [l2j • 

Taking this comparison beyond linear response re- 
quires numerical solution of (|ll[) - (|15p . For simplicity 
we perform calculations for a one-component system of 
hard spheres [T3|. In Fig. [1] we show a xy under shear and 
(cr xx —a yy ) under planar elongation as a function of Peclct 
number, Pen = 7c? 2 / Dq, where 7 is the strain rate, for 
various packing fractions around the glass transition <p c . 
A detailed comparison of theoretical shear flow curves 
with experiment can be found in Q; we do not know of 
similar experimental data on elongational flows but hope 
our work may stimulate future studies. The lower left 
panel shows the Trouton ratio (a xx — o~ yy )/o~xy as a func- 
tion of Pco- Upon entering the nonlinear (shear thinning) 
regime, this ratio is significantly reduced below the linear 
response value, indicating enhanced strain thinning un- 
der elongational flow relative to shear. Moreover, as <p c is 
approached from below, the linear response regime moves 
to lower values of Peo and disappears entirely on cross- 
ing the glass transition. The resulting plateau in the flow 
curve identifies a dynamical yield stress. The Trouton ra- 
tio displays an approximately constant (noninteger) value 
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xy under shear (full lines, 
yy under planar 



FIG. 1: Steady state stress a 
Kij =7 S X iSyj) and stress differe 
elongation (broken lines, Kij = 7 (S x iS x j— 5 y iS y j)) asafunction 
of Peo = 'yd 2 /Do, where d is the sphere diameter. Each curve 
is labelled according to the distance from the glass transition 
<f> - cf> c , (a) -1CT 4 , (b) -1CT 3 and (c) 1CT 4 . The inset shows 
a possible realization of planar elongational flow. The lower 
left panel shows the Trouton ratio for such flow while the 
lower right panel shows the ratio for uniaxial elongation (fty = 



>o"xj ^yi^yj /2" 



-S z iS z j /2)). (Legend (a-c) defined above. 



over the entire plateau: its low Peo limit furnishes a non- 
trivial yield condition for the ratio between normal and 
shear stresses. In the lower right panel of Fig. [T]we show 
in addition the Trouton ratio for uniaxial elongational 
flow. We again obtain a nontrivial yield stress ratio, dis- 
tinct from that found under planar elongation. Our re- 



sults for uniaxial and planar extension, alongside those 
for shear [3] thus predict by first-principles theory three 
parameters of a dynamic yield surface for glasses. The 
limiting stress at zero flow rate for a fixed flow geome- 
try identifies a point on a dynamic yield manifold. In 
a system where (as here) isotropic pressure is irrelevant, 
this manifold resides in the space of deviatoric (i.e. trace- 
less) stress tensors, and in the principal stress frame (with 
stress eigenvalues si :2 ,3) can be represented, for instance, 
as a closed curve, /(si — s 3 ,s 2 — S3) = 0. Our new re- 
sults for uniaxial and planar extension, alongside those 
for shear [7], predict from first principles three parame- 
ters of the dynamic yield manifold for a glass. We plan 
a fuller examination of this manifold in future work. 

In conclusion, we have derived from first principles 
a constitutive equation for dense colloidal suspensions 
subjected to an arbitrary time-dependent (but homoge- 
neous) deformation. Appeal to translational invariancc 
identifies appropriate deformation measures for an ap- 
proximate description of the system in terms of density 
fluctuations. Within our treatment macroscopic stress is 
directly connected with structural distortion at the mi- 
crostructural level. The inclusion of non-exponential slow 
relaxation arising from interactions leads to a unified de- 
scription of fluid and glassy states. For steady shear, pla- 
nar and uniaxial elongational flows we find strong nonlin- 
carities in the Trouton ratio as a function of strain rate. 
For colloidal glasses, we predict nontrivial relations be- 
tween dynamic yield stresses in different flow geometries. 
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